NMR study of the S = 1/2 Heisenberg Ladder Cu 2 (C5Hi 2 N 2 ) 2 Cl4: 
Quantum Phase Transition and Critical Dynamics 
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We present an extensive NMR study of the spin-1/2 antiferromagnetic Heisenberg ladder 
Cu2(C5Hi2N2)2Ci4 in a magnetic field range 4.5 - 16. 7T. By measuring the proton NMR relax- 
ation rate 1/Ti and varying the magnetic field around the critical field H c \ = A/gfiB ~ 7.5T, we 
have studied the transition from a gapped spin liquid ground state to a gapless magnetic regime 
which can be described as a Luttinger liquid. We identify an intermediate regime T > \H — H c i\, 
where the spin dynamics is (possibly) only controlled by the T=0 critical point H c i- 

75.10Jm,75.40Cx,76.60.-k 



One of the most fascinating effect in quantum mag- 
netism is perhaps the possibility to realize a T=0 phase 
transition [EJ from a " quantum disordered" (gapped spin 
liquid) ground state to a Luttinger liquid state in one 
dimension (ID) B or to a Neel ordered state in two di- 
mensions (2D) Hp[. There are well known examples of 
gapped spin liquids in ID: for integer-spin Heisenberg 
antiferromagnetic (HAF) chains ||@] or spin-1/2 HAF 
even- leg ladders 0, quantum fluctuations induce a spin 
gap A between a singlet (S=0) ground state and triplet 
(S=l) excited states. An external magnetic field lifts 
the triplet degeneracy and induces a second order T=0 
phase transition at a critical field H c \ = A when the 
lowest branch of the triplet crosses the ground state. At 
this critical field, dynamical properties are defined by 
universal exponents PHlO|] . Above if cl , the ground state 
is magnetic with an algebraic decay of the correlation 
functions. Close to H c \, the low frequency (u> < T) spin 
dynamics should be governed in an universal way by the 
T=0 critical point H c \. In the vicinity of this point, di- 
vergent quantum correlations are cut-off by thermal fluc- 
tuations at a length £t which become the only relevant 
length scale. This is the quantum critical regime (QCR) 
where the temperature is the only energy scale over the 
large intermediate region J > T > \H — H c \\ j|. The 
true challenge to experiments is to observe in a single 
system, all the sequences of the T — > regimes as well as 
the finite temperature critical regime. 

In a gapped HAF, such quantum phase transition can 
be experimentally studied only if the gap A is comparable 
to accessible magnetic fields. The situation is hopless in 
the S=l/2 spin ladder compound SrCu 2 3 (A > 400K) 
jTlJ. On the other hand, it is now well established that 
the organo-metallic compound Cu2(C5Hi 2 N 2 )2Cl4 ]ll| is 
a unique representative of HAF S=l/2 ladders with small 
exchange constants (J± sw 13. 2K, Jn w 2.5K) and a spin 



gap (A k=Tr ~ Jx~J\\ = 10.5±0.3K) (l|-[T6) which makes 
the entire phase diagram experimentally accessible. 

This Letter describes a proton ( X H) NMR study of 
Cu2(C5Hi 2 N 2 )2Cl4 in magnetic fields ranging from 4.5 
to 16.7 Tesla. The complete phase diagram (including 
the various temperature crossovers) has been observed 
experimentally. We unambiguously identify three dif- 
ferent regimes when T — > 0: 1) a gapped phase for 
H < H c i pa 7.5T, defined by a spin liquid ground state, 
where the energy gap, deduced from nuclear spin-lattice 
relaxation time (Ti) measurements, is linearly reduced 
by the magnetic field; 2) a magnetic phase with a gapless 
ground state for H c i < H < H C 2 rs 13. 2T, character- 
ized by a power-law divergence of 1/Ti consistent with 
an interpretation based on fermions in one dimension. 3) 
a fully polarized gapped phase above H C 2- At interme- 
diate temperatures, we analyze the crossovers between 
these regimes and give for the first time convincing evi- 
dences that a quantum critical regime is observed when 
A>T>|ff-flci|. 

Proton NMR experiments have been performed with 
the magnetic field direction along the b axis (perpen- 
dicular to the [101] ladder direction) of small single 
crystals (typically 100-200^^ each). We used conven- 
tional pulse spin-echo sequences and frequency-shifted 
(summed Fourier transform) processing fl7|| . In a first 
run (H < 8.7T), a set of five crystals with their b axis 
oriented has been used; for larger fields, one single crystal 
was utilised. While spectra displayed minor variations, 
the absolute values of T\ were found reproducible be- 
tween these two runs. A typical proton NMR spectra is 
displayed in Fig. [l]. Note that the large magnetization at 
16.7 Tesla allows to resolve all 24 proton sites. However, 
since the total width of the spectra follows the magnetiza- 
tion M, this discrimination becomes increasingly difficult 
as M — > 0. Arrow in Fig. |] indicates the line (I) where 
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the hyperfine shift K and T-y have been measured (see 
Ref. Q for a discussion of the proton sites assignment). 

The temperature dependence of the magnetic hyper- 
fine shift K ~ /H at various magnetic fields is shown 
in Fig. |l|. The uniform static (q = 0, u) = 0) susceptibility 
Xo = M / H at different fields is also displayed on the same 
scale using the relation K = (Ai / gb^B)xo with the hy- 
perfine coupling Ai rj 2.9kOe and the g-factor gi, = 2.03 
jl4]]. For H < Hci, both K and xo drop exponentially 
due to the effective gap A^ = A — g^sH between the 
singlet and the lowest triplet state. 

As the critical field H c i is crossed, both K and xo show 
a clear persistent downturn at low temperatures due to 
residual short range antiferromagnetic correlations in the 
intermediate phase. When T — > 0, K and xo go to a finite 
value, which rises continuously as the field is increased 
from H c i to H C 2, following the increase of ground state 
magnetization. Above ff C 2 where M( T=0 ) is saturated 
fuj, M/H should decreases as 1/H (H > H c2 , T — > 
0). This is already visible at T « 1.4K, where K at 
16. 7T is smaller than at 15T. The maximum value of K, 
extrapolated to T = corresponds to (S z ) — \. Above 
T J±, we recover a Curie-like tail for all fields, which 
defines the "classical" (decorrelated) regime. 

We now discuss the spin-lattice relaxation rate 1/Ti 
which is expressed in terms of the magnetic structure 
factors ||2Ctl 



Ti~ 2 



^2F a (q) ■ S a (q,u n ) 



q.a 



S a {q,uJ n ) = J exp(iu n t)dt(S a (q,t)S a (-q,0)) , (1) 

where F a (q) are hyperfine form factors of dipolar origin 
^lf and a = z, ± represents the longitudinal and trans- 
verse components, respectively. 

At low fields and low temperatures, T -C H c \ — H, 
(singlet gapped phase), 1/Ti falls off exponentially with 
a characteristic energy gap = A — g^sH [see Fig. 
2(a)]. This effective gap is represented in the lower graph 
of Fig. 2(a). This result can easily be understood: for 
line (I), it has been shown in a previous study |18| that 
two-magnon processes (near k = it, q ~ 0) in the "in- 
trabranch" channel |p2| dominate the nuclear relaxation. 
It follows that 1/Ti, in this gapped phase, is driven by 
longitudinal correlations S zz (q ~ 0, uj n ) ~ exp(— A?-). As 
we approach the critical field, T w H c \ — H , the analy- 
sis breaks down since competing mechanisms like "stag- 
gered" direct, q = it, processes take over and dominate 
the relaxation |^]. At higher temperatures, T > A, 1/Ti 
is constant, as expected in the classical limit pCj], 

In the magnetic phase (H c i < H < H C 2) Jl4| , 1/Ti 
turns upward at low temperature in striking contrast 
with its behavior in the gapped phase [see top parts of 
Fig. 2(a) and 2(b)]. Very close to iJ cl , a divergence is 



readily visible below 2-3K. As the magnetic field is in- 
creased just above H c \ (7.5 < H < 9T), the divergence 
becomes more pronounced and develops at temperatures 
below 5K. At higher fields, the divergent behavior is re- 
placed by a smooth increase as the temperature is lowered 
below 10K and progressively vanishes as we approach the 
upper critical field H C 2. 

To conclude the T = limit, we discuss the "high 
field phase" which appears above H C 2 = J± + 2J\\ [ pT[ . 
It is apparent [bottom of Fig. 2(b)] that the relaxation 
rate decays exponentially at low temperatures with an 
activation energy A up — g^B{H — H c2 ). Since the T=0 
magnetization is saturated above H c2 , the spin system 
is fully polarized H] and all dimers are in the triplet 
state |tt)- A up can therefore be interpreted as the energy 
gap between the "fully polarized" ground state \FGS) = 
IIi||t)i and the lowest excited states generated by a single 
spin flip. The balance between the energy gain J± + 2J\\ 
when the antiferromagnetic couplings are satisfied and 
the Zeeman term g/J-BH determines the value g^isiil — 
H C 2) for the energy gap. 

At this stage, we propose the experimental phase dia- 
gram schematically shown in Fig. |[ The T = phases 
discussed above are: 1) the gapped spin liquid phase; 2) 
a magnetic phase to be characterized (see below) and 3) 
the gapped polarized phase. 

At finite temperatures, these phases are separated by 
crossover lines: (A) and (B) in Fig. || correspond, re- 
spectively, to the onset of the gapped spin-liquid regime 
and of the gapped polarized regime, both characterized 
by an exponential decay of 1/Tl. Line (C) sets the up- 
per boundary of the " magnetic phase" , where a divergent 
behavior of 1/Ti [~ T~ a ] is observed. At higher temper- 
atures (T > A), we have the classical regime defined by 
1/Ti ~ est ||(J and a Curie-like behavior of K ~ 1/T at 
all magnetic fields [see Fig [l| . What can be said about the 
large intermediate region between these crossover lines? 
On one hand, 1/Ti is found to be nearly T-independent 
in the range 6.6 < H < 9T, above a characteristic tem- 
perature "T " ~ \H - H c i\ [see Fig. 2]. On the other 
hand, K shows a clear field dependence below T < A. 
This is no longer a classical regime: we propose that this 
region corresponds to a QCR where 1/Ti is predicted to 
be almost temperature independent H while, at the same 
time, static properties depart from the classical picture. 

In the final part of this Letter, we discuss the nature 
of the magnetic phase, particularly the origin of the low 
temperature divergence of 1/Ti between H c \ and H C 2- 

Since a 3D field-induced ordering above H c i has been 
observed in specific heat measurements below Tn(H) < 
0.8K p3[ , we first consider critical fluctuations as a pos- 
sible origin for the observed divergence. In this scenario, 
one expects a divergence of 1/Ti in a range 5T « T/v 
above T N |24) with a generic behavior Tf 1 - (T-T N )- V 
(0.5 < v < 1). The exponent v must be independent of 
the magnetic field. We cannot fit our data taking into 
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account the known field dependence of Tn(H) |23| with- 
out releasing this constraint. Moreover, 1/Ti starts to 
diverge at T w 5K, which is at least 5T « 5 — 6T/v above 
TV- From these arguments, it is clear that the onset of a 
3D ordering cannot explain alone the behavior of 1/Ti in 
this temperature range (T > 1.3K). Hence, ID quantum 
fluctuations have to be invoked to explain our results [^5| . 

We then propose an analysis in terms of fermions in 
one dimension (i.e a Luttinger liquid in the regime 
T <C H — H c \. The spin-ladder Hamiltonian can be con- 
verted into a ID interacting spinless fermions model using 
a Jordan- Wigner transformation |26|||[l0| . The spectrum 
consists of two bands with energies huik = ±(A + -§%k 2 ). 
In this mapping, the magnetic field H plays the role of 
the chemical potential /i = g/isiH — H c i) of the fermions 
(at H = 0, fjL lies in the middle of the spectrum). Above 
H c i > 0), the density of fermions n increases with fi 
as n cc M ~ ^JJi p7|,p|. Since there is no gap, direct 
nuclear relaxation processes are allowed and one expects 
an enhancement of the relaxation by an amount related 
to the magnetization (M oc n) of the electronic system. 
Within this picture, the staggered part S±(q « n, oj n ) 
leads to 1/T X ~ J 1- Q ( a = 0.5) when H -> H+ 
At higher fields, n increases and interactions give rise to 
non-universal behaviors of the spin correlations, i.e de- 
pendent of the microscopic details. For instance, Ref. 
pp| | predicts that a stays close to 0.5 above H c \ for a 
ladder in contrast with other S=l/2 gapped systems. 

This picture agrees with our results, at least below 9 
Tesla, and becomes poorer above due, possibly, to large 
interactions between excitations @@. The data in Fig. 
2 show that the divergent term (controlled by the ex- 
ponent a), increases from H c \ and is maximum around 
H 8.5 — 9T. When H — * H C 2, the divergence weakens 
and the exponent a cannot be reliably estimated. An im- 
portant point is that the exponent a is related to the ex- 
ponent rj{= 1 — a) controlling the decay of the spatial cor- 
relation ,(SoS r ) ~ (— l) r |r|~''. Our data are roughly con- 
sistent with exact diagonalization calculations of (SoS r ) 
for Haldane chains in the gapless phase ps|] : both rj and 
a are close to 0.5 at H c \ and H C 2- In between, r\ should 
have a minimum value 77 w 0.3 ]2§| ], meaning a maximum 
of a. To summarize, we have shown that the field de- 
pendent 1/Ti divergence is a ID effect and not an onset 
of the 3D ordering which occurs at lower temperature. 
Even though the exponent a cannot be estimated pre- 
cisely, a value of a = 0.5 at H c \ is consistent with our 
data. 

In conclusion, the magnetic field-temperature phase di- 
agram of the spin-1/2 HAF ladder Cus^HiaNa^CLj 
has been completely explored by probing the low fre- 
quency spin dynamics. Three T — > regimes are iden- 
tified. 1) a gapped spin liquid phase for H < H c \; 2) a 
Luttinger liquid phase, well defined for H c \ < H, char- 
acterized by a field-dependent correlation exponent a. 



3) a fully polarized gapped phase H > H C 2- We em- 
phasize that the exponent a derived from our data is 
in qualitative agreement with a Luttinger liquid picture 
very close to H^. Nevertheless, further studies of the 
field dependence of the exponent a are needed. Finally, 
various temperature crossovers are observed and closely 
resembles theoretical predictions for a " quantum critical" 
regime in a range A > T > \H — H c i\. 

We thank R. Chitra, T. Giamarchi, D.K. Morr, D. Poil- 
blanc and S. Sachdev for stimulating discussions. One 
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FIG. 1. Top panel : 1 H NMR spectrum at f = 710.1MHz 
and T=5.1K. All measurements reported here were made on 
the line marked by an arrow. Bottom panel: Magnetic hy- 
perfine shift Ki as a function of H and T. Bold and dotted 
lines represent macroscopic susceptibility data obtained by 
AC-SQUID magnetometer (h=0.1T) and torsional oscillator 
magnetometer, respectively jyj. The normalization factor 
gives an hyperfine coupling Ai ~ 2.9kOe. 



FIG. 2. Panel (a): Temperature dependence of 1/Ti 
through the critical field H&. The lower part covers the sin- 
glet gapped phase (H < H c i) while the upper part is in the 
magnetic phase above H c \ . The dahsed lines are guides to the 
eyes. Lines in the inset are low-T fits to the field dependent 
energy gap Ah — A — g/is-ff with gb = 2.03. Experimen- 
tal values of Ah are shown in the lower graph (□) together 
with those obtained from x(T) and M(H) (•) [|14|. Lines la- 
beled (A) and (C) correspond to the crossover lines of Fig. 
^| Panel (b): Temperature dependence of 1/Ti through the 
critical field H C 2- The upper part is in the magnetic phase as 
H C 2 is approached while the lower part is in the "fully polar- 
ized" gapped phase (H > H C 2). Lines in the inset are low-T 
fits with an energy gap A up = g^iB{H — H C 2) and gb = 2.03. 
Experimental values of A" p are shown in the lower graph (□) 
together with H c2 (•) O. 



FIG. 3. H-T phase diagram which can be tentatively 
drawn from the present experiment. Dashed lines represent 
crossovers between the different regimes discussed in the text: 
line (D) is the QCR-classical crossover, lines (A) and (B) cor- 
respond to T ~ H c i — H and T ~ H C 2 — H , respectively. Line 
(C) is the QCR to Luttinger liquid (LL) regime crossover. 
The black region corresponds to the 3D ordered phase with 
T N (max) w 0.8K & 
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